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H Æ.: 利用 2019 一 2020 4 5—9 H ECMWF (European Center for Medium-Range Weather Forecast) ,GRAPES-GFS 
(中 国 气 象 局 GRAPES 全 球 数 值 预报 业务 系统 , Global/Regional Assimilation and Prediction System-Global Forecast 
System ) 大 尺度 数值 模式 36h 预报 时效 的 24h 降 水 量 预报 和 加 密 降水 观测 资料 ,基于 CRA 空间 检验 技术 识别 ,分离 
青藏 高 原 东 北 侧 暴 十 目标 (ECMWF 模式 18 个 ,GRAPES-GFS 模 式 11 个 ), 定 量 分 析 了 两 种 模式 对 暴雨 预报 的 空间 


误差 ( 落 区 .强度 和 形态 误差 ) 特 征 ,总 结 大 尺度 数值 模式 在 青藏 高 原 东 北 侧 暴 雨 预报 中 的 适用 性 。 结 果 表 明 :(1) 
两 种 模式 的 降水 预报 形态 误差 占 比 均 最 大 。ECMWF 预 报 的 强度 误差 占 比 最 小 , 落 区 误差 次 之 ,而 GRAPES-GFS 的 
落 区 误差 占 比 最 小 ,强度 误差 次 之 ,(2) 两 种 模式 预报 的 暴雨 落 区 位 置 均 较 实况 偏 西 偏 北 。 暴 雨 大 值 中 心 ECMWF 


ii PUES , Mi GRAPES-GFS {Mia PY o (3) 对 暴雨 十 区 面积 的 预报 两 种 模式 均 明 显 偏 小 ,容易 造成 暴雨 漏 报 ,对 最 大 降 


水 量 和 平均 雨 强 的 预报 CRAPES-CFS 模 式 均 偏 弱 超 过 40% , ECMWF 预报 的 平均 雨 强 偏 弱 11.49% ,最 大 降水 量 偏 强 


1.47%。(4) 两 种 模式 对 上 甘肃 陇 东南 地 区 和 陕西 西南 部 的 暴 月 


果 较 差 。 
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降水 预报 是 现代 天 气 预 报 的 重要 业务 之 一 , 数 
值 预报 为 降水 预报 提供 最 有 价值 的 参考 ,但 由 于 模 
式 设计 等 本 身 原因 及 地 形 等 外 部 因素 ,模式 输出 的 
降水 等 预报 产品 在 时 空 分 布 上 均 有 一 定 的 误差 ， 
此 ,通过 检验 和 评估 模式 性 能 为 预报 员 订 正 模式 预 
报 提 供 背 景 误差 信息 中 。 点 对 点 的 常规 天 气 检验 方 
法 (如 TS 评分 和 ETS 评 分 ) 对 降水 落 区 和 时 间 较 为 
敏感 ,存在 “双重 惩罚 "现象 ,同时 也 会 掩盖 一 些 
对 预报 员 订 正 模 式 预 报 有 参考 价值 的 信息 * ,无 法 
对 业务 预报 中 关注 的 降水 强度 、 范 围 大 小 及 十 带 落 
区 等 属性 给 出 定量 检验 结果 ,也 无 法 定量 分 析 预 报 
误差 来 源 ““。 为 了 获取 更 详细 的 模式 预报 偏差 信 
息 ,气象 学 者 相继 开发 了 CRA(Contiguous Rain Ar- 
ea) 、SAL (Structure Amplitude Location) 、MODE 
(Method for Object-based Diagnostic Evaluation ) 等 多 
PF H ERE VE AY 28 a eT PEA, SALZ Ti] 
检验 方法 以 某 一 范围 内 识别 的 降水 为 目标 ,定量 检 
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预报 效果 较 好 ,而 对 陕 北 及 宁夏 等 偏 北 地 区 的 预报 效 


验 降 水 日 标的 强度 (Amplitude,A)、 位 置 (Location， 
L) 415 (Structure, S) IRZ ,该 方法 中 强度 位置 和 
结构 3 个 分 量 是 相对 的 度量 , 且 分 离 这 3 个 分 量 的 
方法 具有 一 定 的 主观 性 ""。MODE 方 法 在 一 定 区 域 
内 识别 分离 实况 和 预报 降水 场 中 的 目标 对 象 , 比 
较 两 个 场 中 目标 对 象 的 空间 属性 ,在 此 基础 上 计算 
目标 对 象 的 相似 度 特征 ,最 终 获 得 实况 和 预报 场 的 
总 体 空间 相似 度 特征 ,但 MODE 方法 会 因 平 滑 半 
径 滤波 冰 值 等 因素 导致 评估 结果 出 现 差异 ,可 能 
会 引起 一 定 评估 误导 问题 。CRA 方 法 将 总 体 误 
差分 解 为 强度 . 落 区 及 形态 误差 ,提供 模式 预报 误 
差 信 息 来 源 ,具有 明确 的 天 气 学 检验 意义 。 

CRA 方 法 在 国内 外 被 广泛 用 于 模式 降水 预报 
和 雷达 邻近 降水 预报 检验 2” 。Sharma 等 利用 空 
间 检 验 技术 评估 了 UM(Met Office Unified Model) 模 
式 的 预报 偏差 ,发 现 模式 对 小 量 级 降水 预报 能 
IR ,对 较 大 量 级 降水 的 预报 能 力 较 差 , 且 随 着 降水 
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量 级 的 增 大 ,位 移 误 差 增 大 。Yu 等 "同样 发 现 随 着 
降水 量 级 增 大 ,模式 预报 能 力 减弱 ,在 强 降水 预报 
中 ,形态 误差 最 大 ,位 移 误差 次 之 ,此 发 现 与 符 娇兰 
等 "利用 CRA 空 间 检验 技术 对 我 国 西南 地 区 东部 
强 降 水 的 ECMWF (European Center for Medium- 
Range Weather Forecast) 预 报 的 检验 结果 相似 。 王 
新 敏 等 “通过 对 2018 年 8 月 3 次 登陆 台风 暴雨 过 程 
的 多 模式 检验 发 现 ,CRA 方 法 能 更 全 面 详 细 的 评估 
模式 的 误差 来 源 , 且 大 部 分 模式 预报 降水 个 体 的 误 
差 主 要 来 源 于 位 移 误差 。 

青藏 高 原 东北 侧 是 东亚 夏季 风 北 边缘 过 渡 
带 , 受 东亚 夏季 风 系 统 、 西 风 带 天 气 系统 和 高 原 
天 气 系统 的 共同 影响 ,是 气候 变化 的 敏感 区 ”” 。 
降水 量 年 内 分 布 极 不 平衡 ” ,雨季 降水 量变 化 主要 
是 由 强 降 水 变化 引起 ” ,发 生 在 7 一 8 月 的 暴雨 占 全 
年 总 次 数 的 63%” ,暴雨 发 生 与 该 地 区 复杂 地 形 或 
局 地 小 气候 密切 相关 ”。 近 些 年 来 随 着 大 尺度 数值 
模式 预报 性 能 的 不 断 提 升 ,模式 对 降水 的 预报 能 
也 有 提高 ,但 对 降水 量 的 预报 能 力 随 着 降水 量 级 的 
增加 而 减弱 ,对 强 降水 的 预报 受 地 形 影 响 明 显 “*"， 
且 模 式 对 降水 强度 和 落 区 的 预报 误差 也 较 大 。 本 
文 利 用 基于 目标 属性 的 CRA 空间 检验 技术 ,开展 
ECMWF 和 中 国 气象 局 GRAPES 全 球 数值 预报 业务 
系统 (Global Forecast System) 两 种 大 尺度 数值 模式 
对 青藏 高 原 东 北 侧 暴 雨 预 报 的 空间 偏差 检验 及 分 
析 ,重点 分 析 模 式 预 报 的 落 区 误差 .强度 误差 以 及 
模式 对 暴雨 的 预报 倾向 ,以 期 为 预报 员 提 供 模 式 预 
报 较 详 细 的 评 佑 结果 ,从 而 有 针对 性 的 订正 模式 预 
报 , 同 时 为 模式 研发 提供 有 价值 的 模式 降水 预报 偏 
差 信息 。 


1 数据 与 方法 


1.1 数据 来 源 与 处 理 

基于 目前 气象 业务 中 的 应 用 情况 及 资料 的 完 
整 性 ,对 ECMWF 和 GRAPES-GFS 模式 20h 起 报 的 
36 了 预报 时 效 的 24h 累积 降水 量 在 青藏 高 原 东 北 侧 
(32°~40°N 100?~111*E) 的 预报 空间 误差 进行 分 
析 。 为 了 减少 模式 版 本 升级 对 预报 性 能 评估 的 影 
响 ,选用 2019 一 2020 年 5 一 9 月 ECMWEF 和 GRAPES- 
GES 两 种 模式 降水 预报 数据 ,降水 观测 资料 采用 兰 
州 中 心气 象 台 提供 的 同时 段 的 自动 站 (基准 站 、 基 
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本 站 和 一 般 站 ) 和 区 域 自动 站 24 h 累积 降水 数据 。 
ECMWF 和 GRAPES-GFS 模 式 的 分 辨 率 分 别 为 0.125°x 
0.125° Fil 0.25°x0.25° ,采用 基于 变 分 技术 的 客观 插 
值 方 法 3 将 GRAPES-GFS 模 式 预报 数据 和 实况 观 
测 站 点 数据 统一 插值 为 0.125°x0.125° 的 网 格 数据 ， 
通过 实例 分 析 表 明 该 方法 具有 较 好 的 插值 结果 ,能 
够 满足 客观 分 析 中 准确 性 和 平滑 性 的 要 求 ”*。 
1.2 研究 方法 

CRA 方 法 是 一 种 基于 目标 的 定量 降水 检验 方 
法 ,将 连续 十 区 (Contiguous Rain Area, CRA ) 作 为 日 
标 进行 检验 ”。 检 验 目标 由 某 一 等 值 线 为 边界 的 单 
个 降水 目标 或 雨 带 而 非 整 个 降水 场 , 因 此 ,能够 帮 
助 预报 员 更 好 的 理解 模式 预报 误差 来 源 ”*” 。 该 方 
法 首先 对 一 定 区 域内 的 降水 场 用 某 一 降水 阔 值 来 
定义 连续 十 区 ,然后 计算 CRA 范围 内 预报 和 实况 的 
降水 质心 面积 .平均 降水 强度 等 统计 量 。 定 义 模 
式 预报 总 误差 为 原始 预报 与 实况 偏差 平方 和 的 平 
均值 ,平移 后 的 误差 为 平移 后 预报 与 实况 偏差 平方 
和 的 平均 值 ,为 了 分 析 模 式 预报 误差 的 来 源 , 将 预 
报 雨 区 按照 预报 与 实况 之 间 均 方 根 误差 最 小 进行 
平移 ,得 到 平移 后 的 预报 十 区 。 由 此 ,可 将 模式 总 
误差 分 解 为 强度 误差 .形态 误差 和 位 移 误差 ,其 计 
算 方法 分 别 为 :强度 误差 为 平移 后 的 模式 平均 降水 
强度 减 去 实况 降水 平均 强度 的 平方 ,位 移 误差 为 总 
误差 减 去 平移 后 的 误差 ,形态 误差 为 平移 后 的 误差 
减 去 强度 误差 ,具体 计算 方法 见 参考 文献 [2,17]。 
1.3 暴雨 天 气 过 程 和 有 效 CRA 

青藏 高 原 东北 侧 发 生 区 域 性 暴雨 过 程 较 少 , 规 
定 该 地 区 24 h 累积 降水 量 达 到 上 暴雨 及 以 上 量 级 降 
水 [>50 mm…(24h) | 的 自动 站 和 区 域 自动 站 超过 30 
站 的 即 为 一 次 区 域 性 暴雨 天 气 过 程 。 按 此 标准 进 
行 筛选 ,2019 一 2020 年 5 一 9 月 该 地 区 共有 37 次 暴 
雨天 气 过 程 ( 表 1)。 基 于 这 些 暴 十 天气 过 程 ,采用 
CRA 检验 技术 识别 降水 目标 。 考 虑 到 青藏 高 原 东 北 
侧 暴 雨 受 地 形 等 因素 影响 暴雨 点 相对 分 散光 ,同时 
大 尺度 数值 模式 对 局 地 暴雨 的 预报 能 力 有 限 “”， 
对 暴雨 连续 十 区 进行 识别 和 分 离 时 >50 mm 以 上 量 
级 的 连续 十 区 范围 不 宜 过 大 也 不 宜 过 小 ,此 处 要 求 
连续 暴雨 区 覆盖 范围 超过 0.5°x0.5° 时 为 一 个 有 效 
CRA。 在 识别 出 的 CRA 范围 内 计算 和 分 析 实 况 观 
测 和 模式 预报 降水 的 暴雨 质心 .面积 .平均 雨 强 .最 
大 降水 量 以 及 总 降水 量 误差 ,并 分 析 平 移 CRA 后 实 


202201.00117v1 


chinaXiv 


66 干 F 


ChinaXiv 合 作 期 刊 


39 卷 


况 和 预报 的 空间 误差 ,重点 分 析 强 度 误差 和 平移 误 
差 。 图 1 给 出 了 2019 年 6 月 4 日 08:00 至 5 H 08:00 
发 生 在 青藏 高 原 东北 侧 的 一 次 区 域 性 暴雨 天 气 的 
实况 .ECMWF 预报 的 24 h 累积 降水 和 识别 的 CRA， 
以 此 说 明 CRA 识 别 、 分 离 和 平移 过 程 。 主 观 上 可 以 
看 出 暴雨 以 上 量 级 雨 带 主要 位 于 甘肃 陇 东 南 地 区 
的 东南 部 ,模式 预报 降水 较 实 况 偏 西 偏 北 ( 图 la)， 
识别 和 分 离 的 观测 场 和 预报 场 的 有 效 CRA 如 图 lb 


所 示 ( 填 值 *3”) ,整体 结果 与 主观 分 析 一 致 。 将 EC- 
MWF 预报 降水 按 均 方 根 误差 最 小 的 标准 向 东 向 南 
平移 后 ,模式 预报 降水 更 接近 实况 观测 (图 1c)。 
ECMWF 和 GRAPES-GFS 两 种 模式 分 别 识别 出 
18 个 和 11 个 有 效 CRA( 表 1) ,部 分 暴雨 天 气 过 程 无 
有 效 CRA ,主要 由 3 种 情况 导致 ,一 是 模式 对 于 暴雨 
漏 报 ,ECMWF 漏 报 了 4 次 , 约 占 总 过 程 数 的 11%， 
GRAPES-GFS 模式 漏 报 了 11 次 , 约 占 30%; 二 是 模 


(a) 预报 和 实况 暴雨 区 (b) 识别 的 CRA O 预报 东 移 后 预报 和 实况 暴雨 区 
102°E 106°E 110°E 102°E 106°E 110°E 102°E 106°E 110°E 
多 T T T T Z Z Cc T T T T T Z Z C T T T T T Z 
3 ẹ 3 ẹ È 5 
g ee E% 5 
Z. cm on on in M 
$ EE EE z 
on foe) loa] isa) foe) 
x Z g Z Z Z 
mn N N N 
102°E 106°E 110E œ  1025E 106°E WOE S s 102°% 106°E 110E ™ 


注 :图 (a) (eo 中 阴影 区 为 实况 ,等 值 线 为 相应 的 ECMWF 预报 ;图 ID) 填 值 "1.2.3.4.5 代表 5 个 不 同 的 CRA ,其 中 "3 "为 有 效 CRA。 


图 1 2019 年 6 月 4 日 08:00 至 5 日 08:00>50 mm 降水 及 识别 的 CRA 分 布 
Fig. 1 24 hours accumulated rainfall(> 50 mm)on 4 June 2019 and CRA identified from its 


表 1 2019 一 2020 年 5 一 9 月 青藏 高 原 东北 侧 暴雨 天 气 过 程 \ 暴 雨 雨 区 面积 及 有 效 CRA 统计 


Tab. 1 Statistics of rainstorm processes on the northeastern side of Tibet Plateau from May to September in 2019-2020, 


as well as the effective CRA and rainfall area 


pea 暴雨 区 面积 km ee yates 芭 雨 区 面积 km Mesa eee 
2019-05-05 1709 x x 2020-06-16 61783 v v 
2019-05-06 1832 x x 2020-06-25 20154 v x 
2019-06-04 13312 v v 2020-07-10 11445 v v 
2019-06-05 9688 x x 2020-07-11 2786 x x 
2019-06-20 15609 v v 2020-07-23 16981 x x 
2019-06-27 33730 x x 2020-07-24 11289 v x 
2019-07-21 21002 v x 2020-08-03 3623 x x 
2019-07-28 6618 x x 2020-08-04 39610 v x 
2019-08-02 25670 v v 2020-08-05 2159 x x 
2019-08-03 26854 v v 2020-08-06 12724 v x 
2019-08-08 4682 x x 2020-08-07 2598 x x 
2019-08-19 1235 x x 2020-08-10 6267 x x 
2019-08-26 11970 v x 2020-08-12 14230 v v 
2019-09-09 38960 v v 2020-08-14 2200 x x 
2019-09-12 19184 v x 2020-08-15 44730 v v 
2019-09-13 10496 x x 2020-08-16 19860 v v 
2019-09-14 55048 v v 2020-08-17 4709 x x 
2020-05-07 5017 x x 2020-08-23 7802 x x 
2020-06-15 3405 x x 


EV 表示 暴雨 过 程 中 识别 出 了 有 效 CRA ,x 表示 没有 识别 出 有 效 CRA。 
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式 预 报 面积 小 ,例如 2020 年 7 月 23 日 西 太平 洋 副 热 
带 高 压 西 北边 缘 偏 南 风 气流 中 的 暴雨 天 气 ( 图 略 )， 
ECMWF 预报 的 暴雨 点 相对 分 散 且 暴雨 十 区 面积 ) 
重 偏 小 , 落 区 位 置 偏 西北 , 故 无 法 识别 出 有 效 CRA; 
三 是 位 置 偏差 太 大 ,例如 ,2019 年 6 月 27 日 的 暴雨 
天 气 ,实况 暴雨 区 面积 33730 km?, 模 式 预报 面积 
14845 km ,但 由 于 预报 的 暴雨 十 区 较 实 况 空间 位 置 
偏差 太 大 ,因此 ,也 无 法 识别 有 效 CRA。 两 种 模式 
分 别 有 15 次 过 程 因 暴雨 十 区 面积 预报 太 小 或 位 置 
偏差 太 大 而 没有 识别 出 有 效 CRA。 


2 结果 与 分 析 


2.1 青藏 高 原 东 北 侧 暴 雨天 气 特 征 

青藏 高 原 东 北 侧 2019 一 2020 年 5 一 9 月 均 有 区 
域 性 暴雨 天 气 发 生 ( 图 2) ,其 中 7 月 下 名、 8 月 上 名 
和 中 旬 发 生 次 数 最 多 。37 次 暴雨 过 程 雨 区 面积 (分 
散 暴 雨 过 程 选取 连续 雨 区 最 大 区 域 计算 面积 ) 从 
1709 km 增加 到 61783 km?( $2 1) ,变化 幅度 大 ,最 大 
十 区 面积 与 最 小 雨 区 面积 相差 近 60000 km’, EW 
别 出 有 效 CRA 的 暴雨 过 程 中 实况 十 区 均 超 过 
11000 km? ,集中 在 8 月 上 旬 和 中 旬 。5 月 上 旬 的 暴 
雨天 气 过 程 次 数 也 较 多 ,该 时 段 暴雨 局 地 性 较 强 ， 
暴雨 面积 往往 小 于 6000 km? (K 1) ,无 法 识别 出 有 
效 CRA。ECMWF 和 GRAPES-GFS 两 种 模式 因 漏 报 
或 预报 偏差 大 而 无 法 有 效 预报 暴雨 的 比例 分 别 高 
达 51% 和 70% ,二 者 对 青藏 高 原 东 北 侧 的 暴雨 预报 
能 力 整体 不 高 。 

已 有 研究 表明 , 西 太平 详 副 热带 高 压 对 西北 地 
区 暴雨 有 重要 影响 2 。 本 文 37 次 暴雨 过 程 中 ， 
70% 为 西 太平 洋 副热带 高 压 西 伸 北 抬 输 送 的 暖 湿 空 
气 ,与 北方 冷 空气 交汇 ,配合 低层 动力 抬升 作用 产 
生 , 这 种 天 气 形势 产生 的 暴雨 主要 发 生 在 7 一 8 月 。 
另外 约 30% 的 暴雨 过 程 是 高 空 低 槽 配合 低层 动力 
抬升 作用 产生 ,这 种 天 气 形势 的 暴雨 主要 发 生 在 西 
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图 2 青藏 高 原 东北 侧 暴 雨天 气 过 程 逐 名 分 布 
Fig. 2 The dekad distribution of rainstorm processes on the 


northeastern side of Tibet Plateau 


太平 洋 副热带 高 压 远离 亚 欧 大 陆 期 间 。 
2.2 青藏 高 原 东北 侧 暴 雨 预 报 总 体 误差 

ECMWF 和 GRAPES-CFS 模 式 预 报 的 青藏 高 原 
东北 侧 暴 雨 落 区 、 形 态 和 强度 3 类 误差 如 表 2 所 示 ， 
两 种 模式 的 预报 形态 误差 占 比 最 大 ,分 别 为 52.36% 
和 52.58% , ECMWF 的 强度 误差 占 比 最 小 ,为 
20.73%, 落 区 误差 次 之 , 占 26.93%, 而 GRAPES-GFS 
模式 的 落 区 误差 占 比 最 小 ,为 16.33% ,强度 误差 
31.23%。 两 种 模式 的 落 区 误差 占 比 和 强度 误差 占 
比 主 要 分 布 在 10%~50% , ECMWF 预报 的 50% 暴 雨 
过 程 形态 误差 占 比 超过 50%(40% 暴 雨 过 程 形态 误 
差 占 比 超过 25%) , GRAPES-GES 模式 预报 的 形态 误 
差 占 比 相对 较 分 散 , 超 过 30% 的 暴雨 过 程 形态 误差 
占 比 超过 50%( 图 3)。 

较 大 的 形态 误差 可 能 与 青藏 高 原 东 北 侧 复 困 
地 形 有 关 。 从 2019 一 2020 年 5 一 9 月 的 观测 暴雨 点 
分 布 可 以 发 现 (图 4a) ,暴雨 分 布 与 地 形 有 密切 关 
系 。 秦 岭南 茵 .六 盘山 一 关山 南 厚 、 关 中 平原 等 地 
暴雨 发 生 频 次 高 ,另外 , 受 河谷 .喇叭 口 . 小 山体 等 
小 地 形 影 响 的 地 区 暴雨 发 生 频次 也 较 高 ,如 榆林 大 
河谷 、 积 石山 一 太子 山 等 地 区 (图 4a)。 这 与 黄玉 起 


表 2 暴雨 落 区 ,强度 平均 误差 及 落 区 ,强度 ,形态 误差 平均 占 比 统计 
Tab. 2 The average location and intensity error of rainstorm [> 50 mm: (24h)`] as well as the percentage of location, 
intensity and pattern error in the total error 
平均 面积 平均 雨 强 最 大 降水 总 降水 量 落 区 误差 强度 误差 形态 误差 平均 纬 向 ”平均 经 向 平均 纬 向 平均 经 向 
误差 /% ”误差 /% ” 量 误差 /% ”误差 /% HE% KE% 占 比 /% 位 移 误差 /(°) 位 移 误差 /(*) 质心 误差 /(°) 质 心 误差 /(*) 
ECMWF -16.20 -11.49 1.47 3.51 26.93 20.73 52.36 -0.34 0.08 —0.36 -0.11 
GRAPES-GFS -68.33 -43.40 —49 33 -79.72 16.33 31.23 52.58 -0.15 0.10 —0.22 0.03 
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i ear = 东北 侧 暴 雨 雨 区 的 质心 较 实况 均 明 显 偏 西 (图 5a， 
go | = ECMWF Ja 图 5c),ECMWF 预 报 偏 西 0.36°, 偏 南 0.11°,GRAPES- 
GFS 模 式 预报 的 质心 偏 西 0.22° ,但 预报 偏 南 偏 北 概 
3 oF 9  ” 率 相当 ( 表 2) ,表明 两 种 模式 预报 的 暴雨 大 值 中 心 
A aol 1 40 存在 系统 偏 西 现象 ,ECMWF 的 预报 系统 性 偏 南 现 
象 较 明显 ,而 CRAPES-GFS 模式 预报 的 暴雨 大 值 中 

20r J20 心经 向 偏差 不 明显 。 预 报 的 暴雨 落 区 较 实况 的 位 
移 误差 分 布 表明 , ECMWF 预报 大 部 分 个 例 的 落 区 


强度 误差 形态 误差 


落 区 误差 


图 3 青藏 高 原 东 北 侧 暴雨 天 气 过 程 落 区 强度 和 
形态 误差 占 比 
Fig. 3 Box plots of the location, intensity and pattern 
error ration in the total error on the northeastern side 


of Tibet Plateau 


等 中 赵 海 英 等 ”关于 地 形 对 暴雨 影响 的 研究 结 
一 致 。 分 析 逐 个 天 气 过 程 也 发 现 地 形 对 暴雨 预报 
有 重要 影响 ,例如 2020 年 8 月 16 一 17 日 暴雨 天 气 ， 
实况 暴雨 点 主要 分 布 在 六 盘山 一 关山 南部 ,而 模式 
预报 的 暴雨 落 区 在 山地 南部 和 北部 均 有 分 布 (图 
4b) ,没有 很 好 的 反映 出 山地 地 形 迎 风 坡 和 背风 坡 
对 降水 的 不 同 影响 。 此 外 ,榆林 大 河谷 附近 降水 受 
局 地 小 地 形 增幅 影响 达到 暴雨 ,而 模式 预报 为 漏 报 。 
2.3 青藏 高 原 东北 侧 暴 雨 预报 空间 误差 

ECMWF Ail GRAPES-GES 模式 预报 的 青藏 高 原 


较 实况 偏 西 偏 北 (图 5b) ,平均 偏 西 0.34° ,平均 偏 北 
0.08"( 表 2) ,大 部 分 个 例 偏 北 的 误差 在 1° 以 内 ,误差 
分 布 比较 集中 (图 $b)。 对 所 有 暴雨 天 气 过 程 的 环 
流 形势 进行 分 析 ,发 现 位 移 误 差 的 成 因 主要 有 两 方 
面 , 一 是 模式 预报 的 高 空 系统 相对 实况 偏 慢 , 二 是 
实况 暴雨 多 出 现在 切 变 线 右 侧 的 南 风 中 , 而 模式 预 
报 多 集中 在 切 变 线 或 低 涡 附近 动力 作用 明显 的 区 
域 。 例 如 2019 年 6 月 4 一 5 日 (图 6a) 和 2019 年 8 月 
3 一 4 日 (图 6b) 的 暴雨 天 气 ,预报 暴雨 区 相对 于 实况 
明显 偏 西 偏 北 ,前 者 模式 预报 高 空 系 统 相 对 实况 偏 
慢 , 暴 雨 区 主要 在 700 hPa 低 涡 切 变 线 附近 ,而 实况 
暴雨 出 现在 700 hPa 切 变 线 右 侧 南 风 中 ,后 者 模式 
预报 的 700 hPa 切 变 线 .水汽 辐 合 区 较 实 况 偏 西 偏 
北 , 对 应 的 暴雨 落 区 同样 较 实况 偏 西北 。GRAPES- 
GFS 模 式 预报 的 落 区 同样 较 实 况 偏 西 偏 北 (图 5d)， 
平均 偏 西 不 足 ECMWF 的 一 半 ,平均 偏 北 位 移 与 EC- 
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图 4 青藏 高 原 东 北 侧 2019 一 2020 年 5 一 9 月 暴利 


注 : 图 4b) 为 2020 年 8 月 16 日 08:00 至 17 日 08:00 暴 雨 .16 日 20:00 风 场 (m:s"') 和 比 湿 (g+kg')。 


次数 (a) 及 2020 年 8 月 16 日 暴雨 实况 .预报 (b) 


Fig. 4 Number of rainstorm distributed from May to September of 2019-2020 on the northeastern side of Tibet Plateau (a), 


observation on 16 August 2020 and its corresponding ECMWF forecasts (b) 
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图 5 ECMWF FU GRAPES-GFS 模式 暴雨 预报 较 实况 十 区 的 质心 位 置 误差 (ac) 和 位 移 误差 (b .d) 散 点 分 布 
Fig. 5 Mass center location bias (a, c), and displacement errors (b, d) between observation and forecasts 
and a, b are for ECMWF, c, d are for GRAPES-GFS 
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注 : 图 四 为 2019 年 6 月 4 日 08:00 至 5 日 08:00 暴 雨 .4 日 20:00 700 hPa 风 场 (ms) 和 比 湿 (g:kg); 图 (4b) 为 2019 年 8 月 3 日 08:00 至 4 日 08:00 暴 
雨 .3 日 20:00 700 hPa 风 场 (ms 和 比 湿 (g:kg)。 

图 6 2019 年 6 月 4 日 (a) 和 2019 年 8 月 3 日 (b) 暴 雨 过 程 实况 与 模式 预报 
Fig.6 The observation and its corresponding ECMWF forecasts for (a) 4 June 2019 and (b) 3 August 2019 


MWF 相当 。 综 上 ,ECMWF 和 GRAPES-GFS 模 式 预 面积 均 较 实况 十 区 面积 偏 小 ,分 别 偏 小 16.20% 和 
报 的 落 区 较 实 况 均 偏 西 偏 北 ,整体 而 言 ,ECMWF 预 68.339 (#22), H ECMWF 预报 的 十 区 面积 离散 度 较 


报 的 落 区 偏 西 更 明显 。 大 ,而 GRAPES-GFS 模 式 预 报 的 离散 度 很 小 (图 
2.4 青藏 高 原 东北 侧 暴 雨 预报 面积 误差 7a) ,表明 GRAPES-GFS 模式 更 易 漏 报 暴雨 过 程 。 实 


ECMWF FU GRAPES-GFS 模式 预报 的 暴雨 十 区 况 雨 区 格 点 数 小 于 70 个 (面积 约 10937.5 km’) , EC- 
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MWF 预 报 普 遍 较 实况 雨 区 面积 偏 小 (图 8a) ,无 法 识 
别 出 有 效 CRA( 表 1), 格 点 数 大 于 70 个 ,ECMWF 预 
报 的 雨 区 面积 较 实 况 偏 大 或 偏 小 概率 相当 ,此 时 大 
部 分 过 程 能 够 识别 出 有 效 CRA( 表 1), 由 此 说 明 实 
况 暴 雨 雨 区 较 大 时 ,ECMWF 预报 能 力 较 强 ,对 暴雨 
十 区 面积 的 预报 效果 较 好 ,而 CRAPES-CFS 模 式 预 
报 的 暴雨 区 面积 普遍 较 实 况 偏 小 (图 8b)。 
25 青藏 高 原 东北 侧 暴 十 预报 强度 误差 

ECMWF 预报 的 平均 雨 强 较 实况 偏 弱 11.49% 
( 表 2), 其 中 56% 的 CRA 的 平均 雨 强 较 实况 偏 弱 , 尤 
其 当 实 况 雨 强 超过 60 mm 时 偏 弱 更 明显 (图 8c)， 
GRAPES-GFS 模式 预报 的 平均 雨 强 仅 个 别 CRA 与 
实况 相当 ,其 余 均 偏 弱 ,平均 偏 弱 43.40%( 表 2), 当 
平均 雨 强 >60 mm 时 所 有 的 CRA 均 偏 弱 , 且 平均 偏 


弱 超 过 45%( 图 8d)。 对 于 平均 雨 强 ,ECMWF 预报 
的 离散 度 较 大 ,而 GRAPES-GFS 模 式 预报 的 离散 度 
小 (图 7b) ,更易 漏 报 暴雨 。 对 于 最 大 降水 量 而 言 ， 
ECMWF 预报 的 大 部 分 CRA 较 实 况 偏 小 (图 8e) ,但 
因 部 分 CRA 的 最 大 降水 量 较 实况 偏 大 的 幅度 较 大 ， 
从 而 使 预报 的 平均 最 大 降水 量 较 实 况 偏 强 1.47% 
(#2), Mi GRAPES-GFS 模式 对 最 大 降水 量 的 预报 
离散 度 很 小 (图 7c) ,最 大 降水 量 平均 偏 弱 49.33% 
( 表 2) , 且 呈 现 一 致 性 较 稳 定 的 偏 弱 现象 (图 8f) ， 
JE , GRAPES-GFS 模式 对 暴雨 以 上 量 级 降水 的 预报 
能 力 差 ,很 容易 出 现 漏 报 。 

ECMWF 预报 的 总 降水 量 ( 平 均 降 水 乘 以 格 点 
数 )58% 的 CRA 偏 弱 ,但 偏 弱 幅度 较 偏 强 幅 度 小 (图 
8g) ,导致 预报 平均 较 实 况 偏 强 3.51%。 对 于 < 50x 
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图 7 暴雨 的 实况 和 模式 预报 的 箱 线 图 
Fig. 7 The observation rainfall and forecasts characteristics on the northeastern side of Tibet Plateau 
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图 8 ECMWF(a.b.c.d) 和 GRAPES-GFS(e fgh) 模 式 预报 和 实况 暴雨 面积 .强度 散 点 图 
Fig. 8 Scatter of intensity and area for observed and forecasted, and a, b, c, d are for ECMWF, e, f, g, h are 
for GRAPES-GFS model 
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10 mm*m” 的 总 降水 量 ,ECMWF 预 报 较 实 况 偏 小 的 
概率 较 大 ,但 对 总 降水 量 >50x10: mm- mf CRA, 
ECMWF 预报 偏 强 偏 弱 的 概率 相当 。ECMWF 由 于 
对 尺度 较 小 的 暴雨 雨 区 面积 预报 偏 小 ,导致 对 < 
50x10: mm*m ?总 降水 量 预报 偏 弱 。 与 ECMWF 相 
反 ,GRAPES-GFS 模式 预报 的 总 降水 量 一 致 性 偏 小 
且 偏 小 达 79.72%( 图 8h) ,说 明 模 式 内 部 大 气 水 汽 循 
环 整体 偏 干 。 通 过 对 比 暴 十 以 上 量 级 的 各 类 误差 
RIL, ECMWF 预报 较 GRAPES-GFS 模式 预报 与 实 
况 更 接近 。 
2.6 模式 对 不 同 区 域 的 预报 效果 分 析 

H T A ECMWF 和 GRAPES-CFS 两 种 模式 对 
青藏 高 原 东 北 侧 不 同 区 域 的 预报 效果 ,统计 分 析 有 
效 CRA 的 空间 分 布 。 ECMWF 预报 的 暴雨 能 识别 的 
有 效 CRA 主要 分 布 在 甘肃 河东 东部 及 陕西 大 部 (图 
9) ,其 中 , 陇 东南 地 区 及 陕西 西南 部 地 区 CRA 出 现 
频次 最 高 ,其 次 是 陕西 中 部 .定西 .天 水 东部 及 平 庆 
地 区 , 陕 北 .宁夏 白银 南部 及 定西 北部 CRA 出 现 频 
次 低 。18 次 暴雨 过 程 中 有 效 CRA 平 均 位 于 陇 东 南 
地 区 和 陕西 西南 部 ,GRAPES-GFS 模式 预报 结果 与 
ECMWF 的 相似 (图 略 ) ,但 其 有 效 CRA 地 理 位 置 相 
对 分 散 。ECMWF 对 暴雨 较 多 的 陇 东南 地 区 陕西 西 
南部 的 预报 效果 较 好 ,而 对 陕 北 地 区 和 宁夏 等 偏 北 
地 区 的 预报 效果 较 差 。 
2.7 模式 预报 暴雨 的 倾向 分 析 

为 了 核实 ECMWF 和 GRAPES-CFS 模 式 对 青藏 
高 原 东北 侧 暴雨 事件 的 预报 倾向 ,假设 两 种 模式 预 
报 的 暴雨 十 区 面积 误差 不 超过 20% 为 命中 ,大 于 
20% 为 偏 强 ,小 于 -20% 为 偏 弱 ,平均 雨 强 .最 大 降水 
量 亦 按 此 方法 划分 ,两 种 模式 的 统计 结果 如 表 3 所 
示 。ECMWF 预报 的 平均 雨 强 偏 强 和 偏 弱 的 频次 相 
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注 : 虚 线 表 示 有 效 CRA 所 在 的 经 纬度 范围 。 
图 9 ECMWF 的 有 效 CRA 位 置 分 布 

Fig. 9 The location of effective CRA forecasted by ECMWF 


当 ,命中 频次 最 低 ,最 大 降水 量 预 报 偏 弱 的 频次 最 
高 ,命中 和 偏 强 的 频次 相当 ;GRAPES-GFS 模式 预报 
的 平均 雨 强 和 最 大 降水 量 ,分 别 有 3 个 和 2 个 命中 ， 
其 余 均 偏 弱 。 对 面积 的 预报 ,ECMWF 预报 仅 33% 
命中 ,17% 偏 强 , 偏 弱 的 概率 达 50% ,而 GRAPES- 
GES 模式 预报 偏 弱 的 概率 高 达 82%。 因 此 ,两 种 模 
式 均 倾向 于 将 降水 面积 .平均 雨 强 和 最 大 降水 量 报 
弱 ,GRAPES-GFS 模 式 更 容易 出 现 漏 报 。 


3 结论 


利用 2019 一 2020 年 5 一 9 月 ECMWF 和 GRAPES- 
GFS 大 尺度 数值 模式 20h 起 报 的 36h 预 报时 效 的 
24 了 hd 降水 量 预报 和 加 密实 况 降水 观测 资料 ,基于 目 
标 属性 的 CRA 空间 检验 技术 ,分 析 两 种 模式 对 青藏 
高 原 东北 侧 暴 雨 预 报 的 空间 误差 (位 移 误 差 ,强度 
误差 和 形态 误差 ) 和 对 暴雨 的 预报 倾向 ,得 到 如 下 


结论 : 


表 3 暴雨 雨 区 面积 ,平均 雨 强 ,最 大 降水 量 偏差 频次 分 布 


Tab. 3 Frequency distribution of rain area, average rainfall intensity and maximum precipitation 


deviation of rainstorm of > 50 mm: (24h)" 


平均 雨 强 最 大 降水 量 
偏 弱 命中 偏 强 偏 弱 命中 偏 强 

面积 偏 大 ECMWF 0 0 3 0 2 1 
GRAPES-GFS 0 0 0 0 0 0 

命中 ECMWF 0 1 5 0 2 4 
GRAPES-GFS 0 2 0 2 0 0 

偏 小 ECMWF 7 2 0 8 2 0 
GRAPES-GFS 8 1 0 7 2 0 
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(1) 青藏 高 原 东 北 侧 aN 预报 误差 中 ,ECMWF [J]. 气象 科技 , 2013, 41(4): 696-702. [Cui Fen’ e, Wang Yong, Li 
All GRAPES-GFS 模 式 的 IBA AS 4 误差 4 比 最 K, 分 别 达 tk Huijun. Performance verification of coastal torrential rainfall fore- 
cast with several numerical products[J]. Meteorological Science 
52.36% Fil 52.58% , ECMWF 小 a P 
H 的 强度 Rae Fl 比 最 | ,为 and Technology, 2013, 41(4): 696-702. | 

20.73% , 误差 次 之 ， 而 GRAPES-GFS 模式 的 落 治 [4] Davis C A, Brown B G, Bullock R, et al. The method for object- 
区 误差 | one 为 16.19%, 强度 用 差 次 之 。 暴雨 based diagnostic evaluation (MODE) applied to numerical fore- 
落 KMH 误差 与 大 尺度 模式 预报 系统 快慢 等 因素 casts from the 2005 NSSL/SPC Spring Program[J]. Weather & 

AK, 而 形态 误差 可 能 与 该 地 区 复杂 地 形 有 关 。 Forecasting, 2009, 24(5): 1252-1267. 
[5] HCE AE, F, 邵 玲 玲 , 等 . 强 对 流 天 气 预 报 检验 新 方法 在 上 海 


(2) ECMWF 模式 预报 的 暴雨 落 区 位 置 较 实 况 
偏 西 偏 北 , 平 均 偏 西 0.34°, 偏 北 0.08°,GRAPES-GFS 
平均 偏 西 不 足 ECMWTF 的 一 半 ,平均 偏 北 位 移 与 EC- 
MWF 相 当 , 而 ECMWF 预报 的 暴雨 十 区 质心 均 较 实 
况 明 显 偏 西 偏 南 ,平均 偏 西 0.36° , 偏 南 0.11°， 
GRAPES-GFS 平 均 偏 西 约 0.22°, 偏 南 偏 北 概率 相 
当 。ECMWF 预 报 的 暴雨 中 心 偏 西 偏 南 , 位 置 偏 西 
偏 北 ,而 CRAPES-GFS 模 式 的 暴雨 大 值 中 心 同 样 仿 
西 , 但 偏差 小 于 ECMWF ,经 向 偏差 不 明显 。 

(3) ECMWF 和 GRAPES-GFS 两 种 模式 预报 的 
暴雨 十 区 面积 均 较 实况 偏 小 ,但 ECMWF 对 格 点 数 
大 于 70 个 的 十 区 ,面积 预报 偏 大 偏 小 的 概率 相当 。 
ECMWF 和 GRAPES-GFS 模 式 预 报 的 平均 雨 强 偏 
弱 , 分 别 偏 弱 11.49% 和 43.40% ,平均 雨 强 增强 时 , 偏 
弱 现 象 更 明显 。GRAPES-GFS 模式 对 最 大 降水 量 的 
预报 呈现 一 致 性 的 偏 弱 现象 ,平均 偏 弱 49.33% ,而 
ECMWF 因 对 最 大 降水 量 预报 离散 度 较 大 而 平均 偏 
强 约 1.47%。 

(4) ECMWF 和 GRAPES-GFS 两 种 模式 对 甘肃 
陇 东南 地 区 和 陕西 西南 部 的 暴雨 预 报效 果 较 好 ,而 
在 陕 北 及 宁夏 等 偏 北 地 区 的 预报 效果 较 差 ,两 种 模 
式 对 暴雨 的 预报 倾向 表明 其 对 暴雨 十 区 面积 的 预 
报 明显 偏 小 ,容易 造成 暴雨 漏 报 , 但 ECMWF 对 暴雨 
的 预报 能 力 强 于 GRAPES-GFS 模 式 。 
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Spatial error characteristics of rainstorm forecasts of large-scale numerical 


model over the northeastern side of Tibetan Plateau 
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Abstract: Rainstorm targets were identified from May to September in 2019-2020 over the northeastern side of 
Tibetan Plateau by using the Contiguous Rain Area technique. There were 18 rainstorm targets in European 
Center for Medium-Range Weather Forecast (ECMWF) and 11 in Global/Regional Assimilation and Prediction 
System-Global Forecast System (GRAPES-GFS) mode. The spatial error characteristics—including displacement, 
intensity, and pattern error—of the two models for rainstorm forecasts were analyzed quantitatively, and the 
applicability of the two models in rainstorm forecasts over the northeastern side of Tibetan Plateau is here 
summarized. The results of the study showed that: (1) The pattern errors of rainstorm forecasts in ECMWF and 
GRAPES-GFS accounted for the largest proportion of the total errors. ECMWF had the smallest intensity error 
ratio for rainstorm forecast, followed by displacement errors, while the GRAPES-GFS model had the smallest 
displacement errors, followed by intensity errors. (2) The rainstorm areas in the two models shifted more 
northwestward compared with observations. The forecast of ECMWF for the rainstorm center was 
southwestward, while that of GRAPES-GFS was only westward. (3) The forecast of rainstorm area produced by 
the two models was obviously too limited, meaning that it could easily fail to predict rainstorm events. The 
GRAPE-GFS model underestimated both maximum precipitation and rainfall intensity by more than 40%, while 
ECMWF under- forecasted about 11.49% for the average rainfall intensity and overestimated about 1.47% for 
maximum precipitation. (4) The two models showed a better forecasting ability for rainstorm in the southeastern 
Gansu and southwestern Shaanxi areas, and a poorer one in the northern Shaanxi, Ningxia, and other northern 
areas. 
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